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Abstract 

We consider moments of the non-diagonal vector, axial-vector, scalar and 
pseudo-scalar current correlators involving two different massive quark flavours up 
to three-loop accuracy. Expansions around the limits where one mass is zero and the 
equal-mass case are computed. These results are used to construct approximations 
valid for arbitrary mass values. 
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1 Introduction 



Current correlators are building blocks for a number of physical quantities. Among them 
is the total cross section of hadrons in electron positron annihilation which is obtained 
from the imaginary part of the vector correlator. The axial-vector correlator leads in an 
analogue way to contributions to the Z boson decay rate. The scalar and pseudo-scalar 
correlators govern the decay rate of Higgs bosons with the respective CP property. 

An important role in the context of current correlators is taken over by moments obtained 
from an expansion of the two-point functions for small external momenta. As far as the 
correlators are concerned where the current couples to the same quark flavour they have 
been used to extract precise values for the charm and bottom quark masses by comparing 
the results for the vector correlator up to four-loop order [IHH] with the moments obtained 
from the experimentally measured cross section (see, e.g., Refs. [9HTT]). Also the moments 
for the other correlators [6HH] can be employed once the experimental data is replaced by 
precise lattice simulations as has been done in Refs. [12l[T3] . Besides the mass values also 
the strong coupling constant can be extracted from the comparison of moments of the 
pseudo-scalar currents obtained both in perturbation theory and with the help of lattice 
gauge theory calculations [T2"|[T5]. 

The focus of the paper lies on moments of heavy-light current correlators, i.e., two-point 
functions involving two quark flavours with different masses. Results for the case where 
one of the masses is zero have been obtained in Refs. [T^l[T5] and [MJUT] to two- and 
three-loop order, respectively. These results have been used in [16j[I7] to reconstruct, in 
combination with high-energy expansions and information about the threshold behaviour, 
approximations valid for all external momenta and quark masses. Applications of these 
results are corrections for single-top production or decay rates of charged Higgs bosons in 
extensions of the Standard Model [THJ[T7]. The three- loop corrections have furthermore 
been used to obtain precise results for the B and D meson decay constants, /# and fz>, 
based on sum rules [T8|fT9] . Also the extraction of fs and fo from lattice gauge simulations 
requires perturbative input. As discussed in Ref. [20] it is particularly desirable to have 
moments of the heavy-light currents for general values of the two quark masses, not only 
for the physical masses of the involved quarks in order to perform simulations for a variety 
of different masses, which then can be extrapolated to the mass values of physical interest. 
We therefore extend the known three-loop results for the moments of the vector, axial- 
vector, scalar and pseudo-scalar current correlators where one of the quark masses is zero 
to the case of two different masses. This is done by computing analytic expansions around 
two limiting cases, the one for equal masses, the other for one of the masses equal to zero. 
Subsequently, approximations are constructed which cover the whole range. 

In the next Section we provide details on the calculation, the results for the pseudo- 
scalar correlator are presented in Section |3j We conclude with a summary in Section HI 
Long analytic formulae are relegated to the Appendix where also numerical results for 
the vector, axial-vector and scalar correlator can be found. All results for the moments 
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can be downloaded from [21] . 



2 Calculation 

Let us in a first step present our notation. We define the vector (v), axial- vector (a), 
scalar (s) and pseudo-scalar (p) currents via 

Ol = ^17^2, 
f = $iifa, 

f = ^ 7 ^ 2 , (1) 

where ipi and ip2 denote the two quark flavours. Using these currents the vector and 
axial-vector correlator is defined through (5 = v, a) 

HV + ») IlVHg^niO? 2 ) = ijdxe^(0\Tjl(x)j s J(0)\0), (2) 

where Il 5 (g 2 ) and U s L (q 2 ) are the transverse and longitudinal contributions, respectively. 
The scalar and pseudo-scalar polarization function reads (5 = s,p) 

q 2 U\q 2 ) = ijdxe^(0\Tj 5 (x)f(0)\0). (3) 

Throughout this paper we consider anti-commuting 75 which is justified as for ipi 7^ ip 2 
only non- singlet diagrams contribute. As a consequence the axial- vector (pseudo-scalar) 
correlator coincides with the vector (scalar) one if either mi or 7712 is zero. Sample 
Feynman diagrams occurring at one-, two- and three-loop order are shown in Fig. [Q 

It is convenient to introduce the dimensionless variables 

* = 4' ( 4 ) 
m( 

x = — , 5 
mi 

which enables us to cast II' 5 (g 2 ) for q 2 — > in the form 

n V) = T^E^- W 

n>-l 

where the bar indicates that the overall renormalization of the polarization function is 
performed in the MS scheme. Note that we have the symmetry relations C„(x) = C^(—x) 
and C v n (x) = C%(-x) (see, e.g., Ref. [22]). 
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Figure 1: Sample diagrams contributing to IT* 5 at one, two and three loops. The thick 
(lower) and thin (upper) lines correspond to quarks with mass m\ and rn,2, respectively, 
and the curly lines represent gluons. 

We perform an explicit calculation only for the transverse part of the vector and axial- 
vector current and reconstruct the coefficients for the longitudinal part with the help of 
(n = -l,...,3) 

Cl, n = (l + xfC p n+1 . (7) 

For n = —1, 0, 1 and 2 and six expansion terms in x and (1 — x) we have checked that 
these equations hold. 

The expansion of the coefficients in a s is given by 




where the arguments x and /i are suppressed. Often it is advantageous to require a 
QED-like normalization IT 5 (0) = where the relation to TL s (q 2 ) is given by 

For the practical calculation we use a well tested set-up which is highly automated in order 
to avoid errors. All Feynman diagrams are generated with QGRAF [23] and the various 
diagram topologies are identified and transformed to FORM [23] with the help of q2e and 
exp [25],[26]. exp also applies the asymptotic expansion for x on a diagrammatic 
level. The expansion around x = 1 leads to a naive Taylor series. In a next step the 
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expressions are passed to FORM where appropriate projectors are applied, the expansions 
in small quantities are performed and traces are taken. Once the analytical expressions 
for each diagram are reduced to scalar integrals the program MATAD [27] is called in order 
to handle vacuum integrals up to three loops. 

The three-loop calculation requires the inclusion of counterterms up to two-loop order 
(see, e.g., Ref. [28]) which we implement in the MS scheme both for the parameters m 1; 
rri2 and a s and the renormalization of the current itself. The latter is equal to the mass 
renormalization for the scalar and pseudo-scalar case and unity for the vector and axial- 
vector current. At that point all remaining poles are local and can be subtracted in the 
MS scheme in order to arrive at ITV). (See Ref. [22] for details.) 

Following this procedure we were able to evaluate the coefficients C n (x) with n < 4 where 
terms up to order x 8 and (1 — x) 9 are included in the expansions around x = and x — 1, 
respectively. 



3 Results 

In this Section we discuss the quality of our expansions in x. At one and two loops we 
can compare to the exact result which provides confidence concerning the validity of the 
three-loop approximations. In the main text we restrict ourselves to the pseudo-scalar 
case. Results for the remaining three correlators are presented in the Appendix and in 
Ref. |2U. 



3.1 One- and two- loop results for the pseudo-scalar moments 

The exact dependence on x of the one- loop results for the pseudo-scalar moments C^(x) 
with n < 4 is given by 



c [ T(x) = - 4 (i + y (1 - x + x 2 ) + l x - 

l + 4x + x 2 4x 3 (2 + x 



+ x 



C ° ] ' P{X) ~ ' (1 + zf lx (-l + x)(l+V +2i/ " 
g(o) lP(x) = 2(l-x + x 2 )(l + 4x + x 2 ) _ 8x 3 

1 [X) 3 (-1+x) 2 (1 + x) 4 X (-1 + x) 3 (1+x) 5 ' 

£(o), P / _ l + 4x-7x 2 + 40x 3 -7x 4 + 4x 5 + x 6 ^ 4x 3 (2-x + 2x 2 ) 

2 ~ 6 (-1 + x) 4 (1 + x) 6 X (-1 + x) 5 (1 + x) 7 ' 



Cf^ix) 



1 + 5 x - 14 x 2 + 145 x 3 - 94 x 4 + 145 x 5 - 14 x 6 + 5 x 7 + x 8 
15 (-1 + x) 6 (l + x) 8 
^ 8x 3 (1 — x + 3x 2 — x 3 + x 4 ) 
~ X (-1 + x) 7 (1 + x) 9 ' 
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30 (-1 + x) 8 (1+x) 10 
-398 x 6 + 356 x 7 - 23 x 8 + 6 x 9 + x 10 
+ 30 (-1+x) 8 (1+x) 10 

4 x 3 (2 - 3 x + 12 x 2 - 8 x 3 + 12 x 4 - 3 x 5 + 2 x 6 ) 
( — 1 + X) (1+x) 

where = ln(/i 2 /m 2 ) and l x = ln(x 2 )/2. These results have been extracted from Ref. [TB] 
where n p (g 2 ) is given in analytic form up to two loops. We have checked the results of 
Eq. (fTOl) evaluating directly the moments both exactly in x and restricting ourselves to 
expansions in x n and (1 — x) n . Note that the latter method has to be applied at three 
loops. 

In Fig. [2] we discuss (7„ p (x) for n = 1,2,3 and 4 and show the curves including x°, 
x 7 and x 8 for the small-x expansion as dashed lines where longer dashes correspond to 
approximations including higher order power corrections. The dotted lines correspond 
to expansions around the equal-mass case where curves including (1 — x)°, (1 — x) 8 and 
(1 — x) 9 (decreasing distance between dots) are shown. The exact result is shown as 
solid curve. One observes that the approximation including corrections up to order x 8 
coincides with the exact result up to x ~ 0.2, for n — 1 even up to x ~ 0.3. On the other 
hand, the approximation based on the 0((1 — x) 9 ) terms is indistinguishable from the 
solid line for x ;> 0.2 — 0.3. These limits are close to the ones obtained by considering the 
approximations including x 7 and (1 — x) 8 terms only. It is obvious that the small gap in 
between can easily be closed with the help of a polynomial interpolation. 

The exact analytic results for the two-loop moments of the pseudo-scalar current are 
already quite lengthy. Thus we exemplify the structure by showing C[ (x) in the main 
text. The results for Cn' P (x) for n^l can be found in Appendix IA.1I For C^' p (x) we 
have 

-m p . . l + 36x- 22x 2 + 36x 3 + x 4 fT . . 2 . o; , Al 8(l + 3x + x 2 ) 

9 (—1+x) (1+x) 9 (— 1 + x) (1+x) 

16 ^x 2 (1 + 6x + x 2 ) 16/ 2 x 3 (18 + 9x 2 -2x 3 + 3x 4 + x 5 ) , . 

I — ^ i- -. t. -, (11) 

9(-l+x) 3 (l+x)° 9 (-1 + x) 4 (1+x) 6 V ; 

where l\- x i = ln(l — x 2 ) and £2 = 7r 2 /6. Note that in Eq. (fTTj) the coefficient of ln(/i 2 /m 2 ) 
is zero. Using the same notation as in Fig. [2] we show in Fig. [3] the results for cl^'^x). 
The same conclusions as in the one-loop case can be drawn: The expansion for x — > 
provides reliable results for x ^ 0.2 whereas the expansion around x = 1 is trustworthy 
for x 0.2 — 0.4. This is true both for the approximations containing the ninth and the 
ones containing only the eighth order in x and 1 — x, respectively. Thus it is promising to 
proceed at three-loop level, where the computation of the exact dependence on x is quite 
complex, in the following way: using the combination of computer programs described in 
Section [2] it is straightforward to compute expansions in x and (1 — x). Afterwards we 
combine the expansions applying a simple interpolation. 
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Figure 2: One-loop contribution to C v n . The dashed and dotted lines correspond to 
the expansions in x and (1 — x), respectively, where the lowest and the two highest 
approximations are shown. The solid line represents the exact result. Note that for n = 4 
the curve including terms up to order (1 — x) 8 and the one including terms up to order 
(1 — x) 9 lie on top of each other in the shown interval. 



The approximations for moments of the vector, axial-vector and scalar correlators have 
the same quality as in the pseudo-scalar case, i.e. they approximate the exact result in the 
same regions of x. We refrain from an explicit discussion in this paper, however, provide 
the figures corresponding Fig. [2] and [3] in Appendix [B] 



3.2 Three-loop results 

Let us in a first step present the expansions around x = and x = 1. Since they are quite 
lengthy we show in the main part of the paper only results for the first moment adopting 
fi 2 = m\ to exemplify the structure of the result. Furthermore, we restrict ourselves to the 
first three terms in the expansion. The complete expressions for general \x and the results 
for the moments with n = —1, 0, 2, 3 and 4 are provided in the Mathematica file [21]. For 
x —7- we have 

, 2) = 35 1055 _ 88 _ 134 2 105 
1 v ' 16 162 s 9 s 27 S 27 4 






Figure 3: Two-loop contribution to C%. The same notation as in Fig. [2] is adopted. 
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Figure 4: Three-loop contribution to The dash-doted line represents the result based 
on the fit as described in the text. For the rest the same notation as in Fig. [2] is adopted. 



where D 3 « -3.0270, S 2 « 0.2604, S e 2 « 7.8517, T[ « -24.2089 [27] and Cs « 1.2021. 
rz/ is the number of massless quarks and the total number of quarks is rif = n\ + 2. For 
ir — > 1 one obtains 



C{ 2) > p (x) = - 



13139 6977 



+ (1 

+ rii 



2592 1728 
s /8267 



Ca+(1 
6977 



/ 14229845 13342139 



V2592 
25 2 



162 81 



(1 



1728 



x) 



289 
4860 



Cs + (1 



93312 
-x? 



:i-xY + 



103680 
72234997 
933120 
409 



4860 



1 



13760759 
207360 ^ 

+ G((l-x) 4 ) 



(13) 



In Fig. H] we show Cn^' p as a function of x using the same notation as at one and two 
loops. We have chosen n\ = 3 which corresponds to a massive charm and bottom quark. 
We again observe that the approximations based on the x 7 terms and the one including 
x 8 terms coincide up to x~ 0.1 — 0.2. Furthermore, also for the expansion around the 
equal-mass case the two highest approximations overlap for x ^ 0.2 — 0.3, depending on n. 
From the experience gained from the one- and two-loop considerations we can conclude 
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that in these regions our expressions constitute perfect approximations to the exact result. 
Furthermore, it is again straightforward to obtain approximations which are valid for all 
x values by a simple interpolation. In order to obtain handy expressions we perform a 
polynomial fit which is valid for the intermediate region x G [0.1,0.5]. Allowing also for 
half-integer exponents of x we obtain 



,appr\ 



r*(2),p ( > 

L/ 2,apprl J/ > 



M2),P ( > 
^.appr^; 



t -^4,appiA x ; 



0.382 + 22.191 x 1/2 + 5.126 x - 138.613 x 3/2 + 149.094 x 2 - 39.042 x 3 



0.284 + 21.369 x 1/2 - 45.897 x + 11.032 x 3/2 + 26.207a; 2 - 13.199 x 3 



-0.307 + 15.358 x 1/2 - 50.366 x + 61.164 x 3/2 - 27.39 x 2 + 1.475 



x 



0.129 + 4.414 x 1 ' 2 - 16.562 x + 20.537 x 3/2 - 8.586 x 2 - 0.031 x 3 . (14) 



These results are valid for n\ = 3 and are shown as dash-dotted lines in Fig. HI We applied 
the same procedure to the one- and two-loop moments. The comparison with the exact 
result provides an estimate for the accuracy of the approximations in Eqs. (I14p to better 
than 5% for x G [0.1,0.5]. For x < 0.1 or x > 0.5 the corresponding expansions should 
be used for the evaluation of the moments since in these regions they provide excellent 
approximations to the (unknown) exact result. Analogue results for the vector, axial- 
vector and scalar correlators can be found in [21] where again rii = 3 has been chosen. If 
necessary it is straightforward to obtain formulae for other values for n\. 



4 Summary 

We have computed three-loop QCD corrections to the moments of the current correlators 
formed by quark fields with different masses. Our final expressions are based on expansions 
around the known limits where either one mass is zero or both masses are equal. We 
present results which are valid for arbitrary quark masses by combining the expansions 
and a simple polynomial fit in the intermediate region. 

In the main part of the paper we concentrate on the discussion of the pseudo-scalar 
correlator. However, numerical results for the vector, axial-vector and scalar correlator 
can be found in the Appendix and analytical expressions are provided in a Mathematica 
file which comes together with this paper. 

One- and two-loop results with exact quark mass dependence and analytic results for 
the expansions at three-loop order are presented in Mathematica format on the internet 
page [2T] . Also the exact renormalization scale dependence of the three- loop moments are 
provided and numerical approximations valid for arbitrary m 2 /mi. 

The results obtained in this paper constitute important input for lattice calculations in 
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the context of semileptonic and leptonic B meson decays. For such simulations non- 
perturbative Z factors for heavy-light currents are needed which can be extracted from 
the corresponding current-current correlators. Whereas in the case of the B s meson 
expansions for small strange quark mass might be sufficient this is not true in the case 
of B c since m c /mb « 0.23 is in a region where the small-mass expansion already breaks 
down. 
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A Results for the moments of the pseudo-scalar cor- 
relator 

In this Appendix we want to complete the results for the moments of the pseudo-scalar 
correlator. The one-loop expressions can be found in Eq. (HI) and x) is given in 

Eq. dHJ). In Subsection I A. II we provide the results for Cn' p (x) with n = —1,0,2,3 
and 4. As far as the three-loop results are concerned we list in Subsection IA.2I the exact 
x-dependence of the /i-dependent terms which complement the results in Eq. (I14j) . 



A.l Exact two-loop results 
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A. 2 Analytic three-loop results 

Let us decompose the the three-loop coefficients as 
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with 1^ = ln(/i 2 /mf). For n\ = 3 one can find approximations for Cn^' P {x 

Eq. dH]). The exact x-dependence of C^\' p (x) and C^ 2 2 P (x), which can be obtained using 
renormalization group techniques, reads 
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71 + 81 x - 182 x 2 + 10179 a; 3 - 4314 x 4 + 10179 x 5 - 182 x 6 + 81 x 7 + 71 x a 

1620 (-1 + x) 6 (1 + a;) 8 

2 (1 + 5 x - x 2 + 20 x 3 - x 4 + 5 x 5 + x 6 ) 



[Li 2 (x 2 )+2l x l 1 _ x 2-{ 2 ] 

135 (—1 + x) (1 + x) 

4 x 2 (1 + 413 x - 334 x 2 + 1559 x 3 - 574 a; 4 + 833 x 5 - 44 x 6 + 15 a; 7 + 3 a; 8 ) 



+/ 



135 (-1 + x) 7 (1 + x) 9 
2 4 x 3 (180 - 180 x + 1260 x 2 - 800 x 3 + 1940 x 4 - 535 x 5 + 545 x 6 - 4 x 7 + 5 x 8 + x 9 ) 
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135 (-1 + x) 8 (1 + x) 10 



1 + 5 x - 14 x 2 + 145 x 3 - 94 x 4 + 145 x 5 - 14 x 6 + 5 x 7 + x 8 



45(-l + x) D (l + x) f 



l 8x 3 (1 - x + 3x 2 - x 3 +x 4 ) 
~ X 3 (-1 + x) 7 (1 + x) 9 
59 (1 + 6 x - 23 x 2 + 356 x 3 - 398 x 4 + 956 x 5 - 398 x 6 + 356 x 7 - 23 x 8 + 6 x 9 + x 10 ; 

60 (-1 + x) 8 (1 + x) 10 
, 118x 3 (2-3x + 12x 2 -8x 3 + 12x 4 -3x 5 +2x 6 ) /43 

IX (-l + xfd + x) 11 + {Y + ni 

196 + 609 x - 2127 x 2 + 67404 x 3 - 61321 x 4 + 180630 x 5 



(24) 



(25) 



4860 (-1 + x) 8 (1 + x) 10 



-61321 x 6 + 67404 x 7 - 2127 x 8 + 609 x 9 + 196 x 
4860 (-1 + x) 8 (1 + x) 10 



™10 



[Li 2 (x 2 ) +2l x k_ x 2 -f a ] 



1 + 6x - 2x 2 + 54x 3 - 18x 4 + 54x 5 - 2x 6 + 6x 7 



lx X 



135 (-1 + x)' (1 + x)" 
2 + 2412 x - 3130 x 2 + 16758 x 3 - 11358 x 4 + 22158 x 5 
135 (-1 + x) 9 (1 + x) 11 

-6505 x 6 + 5562 x 7 - 214 x 8 + 54 x 9 + 9 x 
135 (-1 + x) 9 (1+x) 11 



10 



+ 2 ll x 3 



360 - 540 x + 4500 x 2 - 4400 x 3 + 13200 x 4 - 7055 x 5 
135 (-1 + x) 10 (1 + x) 12 



+ - 



43 



9810 x 6 - 2169 x 7 + 1476 x 8 - 5 x 9 + 6 x 10 + x 11 
135 (-1 + x) 10 (1 + x) 12 
2 l x x 3 (2 - 3 x + 12 x 2 - 8 x 3 + 12 x 4 - 3 x 5 + 2 x 6 
(-1 + x) 9 (1+x) 11 



(26) 



1 + 6 x - 23 x 2 + 356 x 3 - 398 x 4 + 956 x 5 - 398 x 6 + 356 x 7 - 23 x 8 + 6 x 9 + x 10 

60 (-1 + x) 8 (1 + x) 10 



(27) 
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Figure 5: One-loop contribution to C^. The same notation as in Fig. [2] is adopted. 

B Results for vector, axial-vector and scalar correla- 
tor 

In this appendix we present the results for the moments C^(x) for n = 1, . . . , 4 of the 
vector, axial-vector and scalar correlator. For the analytic results and numerical approxi- 
mations for the three-loop moments in the intermediate x-region we refer to [2TJ and show 
in the following the corresponding numerical result in analogy to Figs. (2J [3] and HI 
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Figure 7: Three-loop contribution to C%. The same notation as in Fig. @] is adopted. 
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Figure 10: Three-loop contribution to C" The same notation as in Fig. [4] is adopted. 
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Figure 12: Two-loop contribution to C*. The same notation as in Fig. [2] is adopted. 






Figure 13: Three-loop contribution to C*. The same notation as in Fig. 0] is adopted. 
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